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We report quantum Monte Carlo (stochastic series expansion) results for the transition from a 
Mott insulator to a dimerized Peierls insulating state in a half-filled, ID extended Hubbard model 
coupled to optical bond phonons. Using electron-electron (e-e) interaction parameters corresponding 
approximately to polyacetylene, we show that the Mott-Peierls transition occurs at a finite value of 
the electron-phonon (e-ph) coupling. We discuss several different criteria for detecting the transition 
and show that they give consistent results. We calculate the critical e-ph coupling as a function of 
the bare phonon frequency and also investigate the sensitivity of the critical coupling to the strength 
of the e-e interaction. In the limit of strong e-e couplings, we map the model to a spin-Peierls chain 
and compare the phase boundary with previous results for the spin-Peierls transition. We point out 
effects of a nonlinear spin-phonon coupling neglected in the mapping to the spin-Peierls model. 

PACS numbers: 75.40.Gb, 75.40.Mg, 75.10.Jm, 75.30.Ds 



I. INTRODUCTION 

Almost half a century ago, Peierls demonstrated that 
a one-dimensional (ID) metal coupled to an elastic lat- 
tice could exhibit an instability towards a lattice distor- 
tion of wave vector q = 2kp} This leads to a gap in 
the electronic spectrum at the Fermi energy and, for the 
case of a half-filled band, the ground state is dimerized. 
Experimentally, the Peierls instability can be observed 
in a wide range of quasi- ID materials, e.g., conjugated 
polymerS)^ organic charge transfer salts;^ MX salts;^ 
and CuGeOs.^ To explain quantitatively the properties 
of these materials, several different models extending 
beyond the original non-interacting Peierls model with 
a classical lattice have been proposedi^ These include 
the Su-Schrieffer-Heeger (SSH) model,7&2iiMi the Hol- 
stein modelfi2iiii4 various Peierls-Hubbardi^iiSiii and 
extended Peierls-HubbardiSii^iiSiSSiSi models, as well as 

spin-Peierls models,^^.^3,24,25,26,27,28,29,30,31,32 

The Peierls instability is well understood in the static 
lattice limit (adiabatic phonons) in the absence of 
electron-electron (e-e) interactions. The ground state 
is a Peierls state for arbitrarily small electron-phonon 
(e-ph) coupling. However, the quantum lattice and e- 
e interaction effects are still not completely understood. 
For many quasi-lD materials the zero-point fluctuations 
of the phonon field are comparable to the amplitude 
of the Peierls distortion,^ and this has spurred a large 
number of studies of quantum fluctuations in SSHjSiifliii 
and Holstein,^^'^^ models. Recent numerical studies have 
shown that quantum fluctuations destroy the Peierls in- 
stability for small e-ph coupling and/or large phonon 
frequency in both the spinlessS^ and spin-^ Holstein 
models^ at half-filling. For the SSH model, Fradkin 
and HirschiS carried out an extensive study of spin-^ 
(n = 2) and spinless {n — 1) fermions. In the anti- 



adiabatic limit (vanishing ionic mass), they mapped the 
system to an n-component Gross-Neveu model, which 
is known to have long-ranged dimerization for arbitrary 
coupling for n > 2 but not for n = 1. For non-zero ionic 
mass, a renormalization group analysis showed that the 
low-energy behavior of the n — 2 model is still governed 
by the zero mass limit of the theory. This implies that 
the spin-i SSH model (but not the spinless model) has 
a dimerized ground state for arbitrarily weak e-ph cou- 
pling. Early numerical calculations are also consistent 
with this scenario.^" It should be noted, however, that a 
dimerized state was also predicted for any non-zero e-ph 
coupling in the spin-i Holstein models ^'^^ for which more 
recent large-scale calculations have instead indicated a 
non-zero critical coupling^ 

Independent-electron models, such as SSH and Hol- 
stein, are important from a theoretical standpoint but 
are not sufficient to account quantitatively for the exper- 
imentally observed properties of real materials. For that 
e-e interactions have to be included in addition to the 
e-ph couplings^ The interplay among the different inter- 
actions gives rise to a rich variety of broken-symmetry 
ground states as well as low-energy electron-lattice ex- 
citations like solitons, polarons, bipolarons, etci^ At 
half-filling, on-site (Hubbard) interactions open a charge 
gap,^'^ and in the absence of e-ph couplings the system is 
then a Mott insulator with algebraically decaying (l/r as 
a function of distance r) spin-spin correlations. Hence, 
the Peierls transition in this case is accompanied only by 
the opening of a spin gap, the charge gap already gen- 
erated by the e-e interactions. Longer-range e-e inter- 
actions can destroy the Mott state, however, and hence 
must affect also the Mott-Peierls transition. Without 
phonons, even in the simplest half-filled extended Hub- 
bard model with only on-site ([/) and nearest-neighbor 
{V) interactionS)2^i2^ some features of the phase diagram 
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are still controversiali22i22i22i2S Adding e-ph interactions 
further increases the complexity of the problem, and the 
determination of the phase diagram remains a very chal- 
lenging problem. 

Zimanyi et aim^ have investigated models with both e-e 
and e-ph interactions using "g-ology" and RG techniques. 
They showed that the ground state has a spin gap if the 
combined backscattering amplitude gf — gi{ujo)+gi{uJo) 
is negative, where gi{ujQ) is the contribution from e-e in- 
teractions and gi^LUo) < comes from the e-ph interac- 
tions. Thus, in the extended Hubbard model, if the bare 
couphng gi ^ U ~ 2V is positive, then > for small 
values of the e-ph coupling and there is no spin gap. The 
transition to a Peierls state occurs only when the e-ph 
coupling exceeds a critical value and gj becomes nega- 
tive. It should be noted, however, that the conventional 
scenario for the behavior close to the line U = 2V has 
recently been challenged^Si^SiSS and a Peierls-like bond- 
ordered state most likely appears close to = U/2 even 
in the absence of e-ph couplings. For the pure SSH model 
([/ = y = 0), (7i is zero and gf is negative for any non- 
zero e-ph coupling. This implies a Peierls ground state 
for arbitrarily small e-ph coupling, in agreement with the 
earlier results of Fradkin and HirschiiS 

In the limit of strong on-site e-e interactions, which 
inhibit doubly occupied sites, a half- filled system can 
be mapped to a spin-phonon model, which also can 
undergo a dimerization (spin-Peierls) transitioni22i24 
Extensive studies of spin-Peierls models in recent 
yeara^^i^^'^^i^^'^°i'^^i^^ have largely been spurred by the 
discovery of a spin-Peierls transition at unusually high 
temperature (14 K) in CuGeOaA Several different cal- 
culations, for different types of spin-phonon couplings, 
have shown that the transition occurs only above a finite 
spin-phonon coupling in the presence of finite-frequency 
phonons.'^^'^^'^"-'^^''^^ This is in contrast to the adia- 
batic limit, where dimerization occurs for infinitesimal 
couplingi^i 

Numerical studies of models with both e-ph and e-e 
interactions, in which the charge degrees of freedom are 
retained (Peierls-Hubbard and extended Peierls Hubbard 
models), have addressed the effect of interactions on the 
dimerization amplitude-^^- -- and the excited states?2i De- 
tailed studies of the phase diagrams have in the past 
been limited by the small lattice sizes accessible when 
both e-e and e-ph interactions are included:-^ The situa- 
tion is rapidly improving, however, as modern quantum 
Monte Carlo22i2L^ and density matrix renormalization 
groupSi*^ methods can now access models with both c- 
e and e-ph interactions on chains with several hundred 
sites. 

Here we consider a ID extended Hubbard model with 
on-site (U) and nearest-neighbor {V) e-e interactions and 
couple it to dispersionless optical bond phonons via mod- 
ulation of the electron kinetic energy. We study the tran- 
sition from a Mott insulating state, with dominant spin- 
spin correlations, to a Peierls (dimerized) spin-gapped 
state. Since the parameter space of this model is rather 



large, with a bare phonon frequency (wq) 3'iid an e-ph 
coupling (a) in addition to the e-e couplings, we have lim- 
ited our study to a physically reasonable ratio U/V ^ 4 
of the e-e parameters. 

In Sec. II we define the model and the various physi- 
cal quantities that we have calculated using a quantum 
Monte Carlo method (stochastic series expansioniS) . In 
Sec. HI we discuss several signals that we have used to 
detect the Mott-Peierls transition. In Sec. IV we present 
the phase diagram in the (wq, Q;)-plane for a fixed value 
of the on-site interaction U that has previously been used 
in models of polyacetylene. We also discuss the effects of 
varying the e-e interaction strength at fixed coq. We map 
the model to a spin-Peierls model for large U, V and com- 
pare the Mott-Peierls boundary with known spin-Peierls 
results. In Sec. V we summarize our results and discuss 
some future prospects. 

II. MODEL AND OBSERVABLES 

The Hamiltonian is given by 

H = -<^(l + a[aJ+ai])(cJ+i,^c,,^-|-4^Ci+i,<,) 

i i 

-hy^(n, -l)(ni+i-l) + tJo5I°^l"^' 

i i 

where a\{ai) creates(annihilates) a phonon on the bond 
between sites i and i -I- 1, c\ ^{ci^cr) is the spin-cr electron 
creation(annihilation) operator, and Ui = ni_i For 
the half-filled band that we study here, the chemical po- 
tential /i = 0. We set the single-electron hopping t to 
unity. For the e-e interactions, we first take the values 
U = 2.5, V = C//4, which have previously been used as 
values corresponding approximately to what is expected 
in polyacetyleneiSS We will also consider other values of 
U, keeping the ratio fixed at U/V = 4. We study the 
system as a function of the bare phonon frequency coq 
and the e-ph interaction a. 

The dispersionless optical phonons we use are different 
from the bare SSH phonons, which have vanishing energy 
for momentum q 0. However, these acoustic phonons 
decouple from the electronic low-energy states involved 
in the Peierls instability and therefore only the optical 
phonons close to q = tt need to be kept.^-J^ Hence, in 
this regard we expect the optical phonons in ^ to be 
equivalent to fully quantum mechanical SSH phonons. 
In the non-interacting limit {U, V ^ 0), the ground state 
should therefore be a dimerized Peierls insulator for any 
non-zero ai^ 

To obtain numerically exact ground state results we 
have used the stochastic series expansion (SSE) quan- 
tum Monte Carlo method"*" for periodic chains with up to 
N = 256 sites. The SSE method is a finite-temperature 
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technique based on importance sampling of the diagonal 
elements of the Taylor expansion of e^^^ , where (3 is 
the inverse temperature; /3 = t/T. Ground state expec- 
tation values can be obtained using sufficiently large (3, 
and there are then no approximations beyond the sta- 
tistical errors. Typically, f3 = 2N or AN was sufficient 
for the quantities presented here to have converged to 
their ground state values. Using the recently developed 
"operator loop" update 4^ the electronic degrees of free- 
dom are treated in the same manner as in the absence of 
phonons.'^'' The phonons are also treated in the occupa- 
tion number basis directly with the SSE representation^ 
(i.e., slightly different from the interaction picture used 
i^or the phonons in in Refs. l29l and li^ . At the (low) en- 
ergy scales that we are interested in here, the number of 
phonons per bond is small (typically < 10) and there are 
no problems in using a truncated basis (the truncation 
can be arbitrarily large in the SSE). 

The Mott and Peierls phases can be characterized us- 
ing the static spin (S) and bond (B) structure factors and 
susceptibilities, as will be further discussed in Sec. III. 
The structure factors are defined by 

^^('?) = ^E^^'''~'^(^fe^f)' (2) 

k,l 

= ]^Y.^''^'''^{KkKi), (3) 

k,l 

where Kj = X]o-=t i^'^j+i a'^j-x' + ^-C-) is the kinetic en- 
ergy operator on the «th bond. The corresponding static 
susceptibilities are given by 

^^(9) = ^E^''^'"'^ f dT{K^{T)Km)- (5) 

k,l •'^ 

Direct evidence for the presence or absence of the spin 
and charge gaps can also be obtained from spin and 
charge stiffnesses pc and ps, which are defined as the 
second derivative of the internal energy per site with re- 
spect to phase factors multiplying the kinetic energy; 
Pes = d'^E{(j)es)/d(l)l^,^'^^' The SSE estimators for all 
these observables can be found in Ref. [sj. 



III. DETECTING THE MOTT-PEIERLS 
TRANSITION 

For our choice of 1^ = C//4, the ground state in the limit 
of zero e-ph interaction (a = 0) is a Mott insulator with 
no spin gap (but finite charge gap) and is characterized by 
a 1/r decay of the staggered spin-spin correlationsi ?^i?7 
The transition to a dimerized Peierls state is marked by 
the development of a staggered kinetic energy modula- 
tion and will hence be signalled by divergent peaks at 
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FIG. 1: Charge (upper panel) and spin (lower panel) stiffness 
constants vs e-ph coupling for several different system sizes 
at ojo = 1 and U = 2.5. 

g = TT in the bond-order structure factor, Eq. Q, and 
susceptibility, Eq. Q. The dimerization is accompanied 
by the opening of a spin gap, the charge gap remaining 
finite. Hence the q = n peak in the spin structure factor 
and susceptibility, Eqs. jSJ and which diverge in the 
Mott phase, become non-divergent in the Peierls state. 

In the adiabatic limit, the system is dimerized for any 
a > 0. We here present several results showing that 
the Peierls transition occurs at a critical coupling etc > 
when ujQ — I and U — 2.5. The phase diagrams discussed 
in Sec. IV are based on the same signals for the transition 
at other e-e couplings and luq. 

Since the charge gap is finite in both the Mott and 
Peierls states, the charge stiffness pc should vanish in the 
thermodynamic limit for all a. The upper panel of Fig. ^ 
shows Pc as a function of a for several system sizes TV. 
As A'^ grows, Pc indeed rapidly converges to zero for all 
a, in agreement to our expectations. The Mott state 
has no spin gap (finite spin stiffness) whereas the Peierls 
state has a finite spin gap (zero spin stiffness). If the 
Peierls transition occurs at a critical coupling ac > 0, 
it should be of the Kosterliz-Thouless type,?- where the 
spin stiffness changes discontinuously from a finite value 
for a < ac to zero for a > ac- The spin stiffness graphed 
in the lower panel of Fig.^shows a jump developing with 
increasing A'^, indicating a critical coupling ac ~ 0.3. 

The spin stiffness data do not easily yield a more ac- 
curate estimate of the critical coupling. As discussed in 
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FIG. 2: Finite-size scaling of the staggered bond-order (upper 
panel) and spin (lower panel) susceptibility for ULjo/t = l and 
several values of a. The dashed lines show the A'^-independent 
behavior expected att the Mott-Peierls transition. 



the context of the spin-Pcierls models* logarithmic cor- 
rections lead to large finite-size effects for a w ac. A 
more accurate estimate can be obtained from the scal- 
ing behavior of the finite-size staggered bond and spin 
susceptibilities)^ It is known from bosonization studies 
that in the Mott phase the equal-time staggered spin and 
bond correlations both decay with distance r as 1/r, up 
to multiplicative logarithmic corrections At the Mott- 
Peierls phase boundary, the log-corrections can be ex- 
pected to disappear)^ and this can be used as a criterion 
for the phase transition. In the dimerized Peierls phase 
the bond correlation function approaches a constant at 
long distances, whereas the spin correlations decay expo- 
nentially. It is convenient to study the associated static 
susceptibilities defined in Eqs. (0J) and ||SJ), which in a 
critical state scale with one power of N higher than the 
structure factors Eqs. |(2J) and ©. In the Peierls phase 
Xs{'^)/N should converge to and xb{t^)/N should di- 
verge, whereas in the Mott phase xs{t^) /N should diverge 
logarithmically and xb{t^)/N should approach zero log- 
arithmically (the log-corrections for spin and bond cor- 
relations are different) . Fig. |2] shows both quantities 
versus In(iV) for several values of a. The expected be- 
havior is indeed observed, and within statistical errors 
both xs(7r)/A^ and xb{'^)/N are independent of N for 
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FIG. 3: Ss{q)/q vs q for several values of a around Oc for 
A*'=128. The curves for a > 0.30 dip below 1 for small q, 
indicating the presence of a spin-gap. 

the largest chains when a w 0.26. 

Additional confirmation of the critical coupling is ob- 
tained by studying the behavior of Ss{q)/q as g — > 0. It 
has been shown^^*^ that if the ground state is gapless in 
the spin sector, then Ss{q)/q ^ I/tt as g ^ 0, whereas 
if there is a spin gap, Ss{q)/q — > 0. Even a very small 
spin gap can be detected this way, since it is in prac- 
tice sufficient to see that nSs{q)/q decays below 1 for 
small q to conclude that a spin gap must be present. In 
Fig. Owe present results for different values of a. The 
curves for a < 0.26 are above I for all q, and the decay 
towards 1 is very slow. The asymptotic approach to 1 
can be expected to be logarithmic4& On the other hand, 
the a > 0.30 curves drop below 1. From these results 
we estimate ac=0.28±0.01, which is compatible with the 
q = TT quantities in Fig. |5] In general, we have found 
that Ss{q)/q, which indirectly signals the opening of a 
spin gap, is the easiest and most reliable way to detect a 
bond-ordered state (see also Ref . . 



IV. PHASE DIAGRAMS 

The critical coupling depends on the parameters of the 
Hamiltonian, in particular, the bare phonon frequency 
ujQ. Using the above criteria for distinguishing the Mott 
and Peierls phases, we have calculated ac as a function 
of uiQ, keeping U — 2.5, V — 0.625. As shown in Fig. 01 
ckc decreases linearly to zero for small loq. This phase 
diagram is hence consistent with the known ac ^ for 
adiabatic phonons. 

For polyacetylene Fradkin and HirscbiS used rescaled 
phonon parameters, which in our units correspond to 
ujo = 0.067, a = 0.052. This point is indicated in Fig.|3| 
and, in accordance with the strong dimerization of poly- 
acetylene, is well within the Peierls phase. 
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FIG. 4: Phase diagram for C/=2.5, 1^=0.625. The squares 
with error bars show the critical e-ph coupling separating 
the Mott (M) and Peierls (P) insulating phases. The cir- 
cle corresponds to phonon parameters previously used for 
polyacetylene^ 



FIG. 5; Phase diagram for ujq — 1,V — f//4. The squares 
with error bars show the critical e-ph coupling. The solid 
curve shows the form ac ~ U"'^. The dashed curve shows the 
transition for the corresponding spin-Peierls modeli— 



As argued above, the e-ph coupling in the present 
model is similar to that in the SSH model for the purpose 
of studying the dimerization transition. This implies that 
in the limit of {U,V) — > 0, we should be able to repro- 
duce previous SSH results. In particular, according to 
Fradkin and HirschjiS ac should be zero even for finite 
frequency phonons. To verify this, we have studied ac as 
a function of {U,V), keeping a fixed ratio of U/V — 4. 
With this ratio the ground state for all [/ is a Mott in- 
sulator with zero spin gap in the limit of vanishing e-ph 
interaction.'^'' We have studied only a single phonon fre- 
quency, = 1 • The resulting phase diagram is presented 
in Fig. 13 The critical e-ph coupling decreases monoton- 
ically with decreasing (U,V), but the smallness of the 
spin gap as {U, V) ^ makes it hard to obtain reliable 
results below U — 0.4. We can therefore not make a defi- 
nite statement about this limit. Nevertheless, our results 
are consistent with a power-law behavior ac ~ with 
7 w 0.3, but a logarithmic form for U ^ can also not 
be excluded. 

For large U, V, and with the ratio V — C//4, the ex- 
tended Hubbard model can be mapped onto the spin-i 
Heisenberg chain with exchange coupling J — 4t^/(t/ — 
V) (this mapping becomes invalid for V ~ where 
phase transitions to bond-ordered and charge-ordered 
phases occuri^^*^Si2i) . Carrying out this transforma- 
tion for the full electron-phonon model the phonon- 
modulated exchange is 



spin-Peierls model considered in, e.g., Refs. 123 and I 



J{Xi 



At'^{l + ax,y 
U -V 



(6) 



where x, 



Under the assumption that the non- 



linear term ~ (cKa;f) can be neglected, which is not a 
priori clear when a is not small, we obtain exactly the 



SP 



= ^^{J + gxi)Si ■ s 



with 



J 



At' 



U-V 



8at' 
U-V 



(7) 



(8) 



For the model Q, an analytic expression for the crit- 
ical spin-phonon coupling g has been obtained for the 
whole range of bare phonon frequencies ujQ/Jm^ The 
form, Eq. (12) of Ref. |3^ is expected to be exact in 
the anti-adiabatic limit, loq/J —>■ oo, which here corre- 
sponds to U —>■ oo. It is in good agreement with numer- 
ical (SSE) results^? for the spin-Peierls transition even 
for frequencies as low as ujq/J = 0.25. In Fig. [3wc com- 
pare our SSE results for the extended Peierls-Hubbard 
model with the spin-Peierls form for U up to 10. The 
transition curve crosses the Mott-Peierls transition curve 
at [/ « 8, and is not in good agreement away from this 
point. The spin-Peierls critical a increases linearly with 
U as U ^ oo, whereas the Mott-Peierls boundary has a 
slower increase with U. The poor agreement for small U 
is not surprising, as the mapping to the spin chain, with 
coupling J given by (jSJ, can only be expected to be good 
for large U , and the form used for the spin-Peierls tran- 
sition curve is not expected to be quantitatively accurate 
for very small oj/J (corresponding here to small U)m^ 
Considering, however, that the analytical form is accu- 
rate for the spin-Peierls model with ojq/J — 0.25; ^^i'^^ 
and that the nonlinear spin-phonon coupling should be 
negligible for small a, the disagreement for the small-JJ 
region in Fig. [S] must be due to the poor correspondence 
between the full electron model and the spin chain with 
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the lowest-order J in Eq. JSJ. The poor agreement for 
U > 8 indicates that the nonUnear spin-Phonon term 
(axi)^ does becomes important as C/ — > oo. 

One effect of the nonUnear couphng term is to renor- 
mahze J: Writing xf — {xf) + A{x'^), the renormahzed 
J is given by 



eff 



4^2(1 



u-v 



(9) 



i.e., Jeff > J. Evaluation the spin-Peierls transition curve 
using Jeff instead of J clearly would reduce ac for given 
U and bring the result closer to the actual Mott-Peierls 
curve in Fig.O However, we have not evaluated (xf) for 
a quantitative test of this effect. In any case, if {xf) is 
large there is also no reason to expect that the remaining 
nonlinear coupling A{xf) can be neglected if a is not 
small. This issue clearly deserves further study. 



V. SUMMARY 



Our phase diagrams are in agreement with what is gen- 
erally expected, but to our knowledge they have not been 
computed quantitatively before. For large e-e couplings, 
we have pointed out the relevance of an effective non- 
linear spin-phonon coupling in the mapping of the Hub- 
bard model to a spin chain. Because of this, standard 
spin-phonon models, where the nonlinear term is not in- 
cluded, cannot be expected to reproduce fully the phase 
diagrams of electron-phonon models. 

The methods that we have used here should also be 
applicable to systems away from half- filling. We plan 
such calculations aimed at studying the stability of the 
soliton lattice, which is formed in doped systems in the 
adiabatic limit, ^'^"■^^ in the presence of finite-frequency 
phonons. A previous quantum Monte Carlo study has 
indicated that the soliton lattice is stableiiS, However, 
open boundary conditions were used for relatively small 
chains, and therefore the issue of whether this is a stable 
phase on an infinite lattice remains to be clarified. 



In summary, we have studied several aspects of the 
phase diagram of an extended ID Hubbard model cou- 
pled to optical bond phonons. We have demonstrated 
that the stochastic series expansion technique^Si^ can 
be used for large electron-phonon chains (here up to 256 
sites) to compute several different quantities that signal 
the opening of the spin gap at the Mott-Peierls transi- 
tion. The spin gap boundary is also in good agreement 
with direct probes of the bond order (i.e., kinetic-energy 
correlation functions and susceptibilities ai q = tt). 
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